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We compared the frequency shift of the ^^Ne Zeeman resonance induced by polarized Rb vapor to 
the shift induced in the ^He Zeeman resonance. The ^He/Rb shift has recently been measured with 
' high precision permitting the conversion of our differential measurement to an absolute value 

\ for the ^^Ne/Rb shift. We report a value of K21 = 31.8 ± 2.8 for the Rb-^^Ne enhancement factor at 

■ a temperature of 124.8° ± 1.9°C. We propose high-precision differential contact shift measurements, 

0^ ' the absolute accuracy of which would be limited by the error in the '^He contact shift. 
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I. INTRODUCTION 



^ 1' Polarization of noble gases by spin exchange with optically pumped alkali metal vapors is an important 
technique for a variety of scientific and technological applications. Polarized noble gases been used for: precise mea- 
surements of the neutron spin structure function 10]; neutron polarizers and analyzers p ] ; magnetic resonance imaging 
studies of gas diffusion in porous media gain media for Zeeman masers and tests of fundamental 

symmetries p^-p^. 

^ New tests of fundamental symmetries using polarized noble gases have recently been proposed. Kostelecky and 

Q Lane recently reported that a ^^Ne/'^He differential magnetometer could be used to test local Lorentz invariance (LLI) 

10. We have proposed to develop a ^^Ne/'^He dual noble gas maser to carry out such an LLI test ||l7|| . 
^>~j. The spin exchange interaction between alkali and noble gas atoms has been extensively described in the litera- 

• ture p| jT^j20| -p^ . During spin exchange collisions, there is significant overlap of the alkali atom's valence electron 
Oh, wavefunction with the noble gas nucleus which comprises a Fermi contact interaction of the form 

^ ■ Hm ^aK-S (1) 

' where S and K are the electron and nuclear spin angular momentum operators, respectively. This interaction leads to 
00 ] exchange of angular momentum between the alkali valence electron and the noble gas nucleus. The valence electron 
. and nucleus experience during the collision a strong magnetic field from the close proximity of their spin magnetic 
moments. In spite of the small fraction of the time spent in collision, the average magnetic fields experienced by 
both the alkali valence electron and the noble gas nucleus can be affected if either the alkali or noble gas species is 
polarized. The resultant change in the Larmor precession frequency for either species, referred to as the contact shift, 
is expressed in terms of an enhancement factor k. k is defined in terms of the total shift in the Larmor precession 
frequencies in co-located spherical distributions of magnetizations, induced both by the contact interaction and by 
the classical magnetic fields generated by the magnetizations, n is the ratio of that total frequency shift (contact 
plus classical magnetization) to the shift that would be induced by the classical magnetization alone, k can be much 
' greater than unity, increasing with increasing Z of the noble gas nucleus [ pO[ . In the limit where the formation of 
O |. alkali-noble gas van der Waals molecules is a negligible part of the spin exchange interaction (true for the experiment 
L| reported in this paper) , a single enhancement factor describes both the contact shift of the noble gas Zeeman splitting 
, by the alkali atoms, and vice versa. [ p3| . The frequency shift experienced by noble gas atoms in the presence of 

• polarized alkali atoms, including both the contact shift and the shift due to classical magnetization fields, is 

S ■ 8^ 

■ - - ' SuJng = {k-1 + foDp) Ing ' -^gSfJ-B [a] (S^) (2) 

Recall the elementary result that is the magnetic field strength inside a magnetized sphere with magnetic 

dipole moment density M; the field is uniform and parallel to the magnetization. foDF is a dimensionless factor 



that accounts for the effect of classical magnetization fields ("distant dipole fields" |18|), and is unity for a spherical 
distribution of magnetization, [a] is the number per unit volume of alkali atoms, and 7„g is the gyromagnetic ratio 
of the noble gas species. For convenience, we define the component of alkali atom magnetization along the magnetic 
field (z) axis, Mq, as 

Ma = n-ia [a] (S,) (3) 
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where 7a = gsi^-B/^ = 27r-1.40 MHz/gauss is the electron gyromagnetic ratio. The alkah nuclear magnetic moments 
do not contribute to the contact shift, but do affect the distant dipole fields produced by the alkali magnetization 
distribution. We ignored these effects on the distant dipole fields, since they were too small to bear on the results 
reported in this work. 

In precision measurements, it is important to be able to monitor the absolute polarization of the noble gases. 
Romalis and Gates jl^ demonstrated accurate polarimetry of "^He using the Rb Zeeman frequency shift due to Rb- 
^He spin exchange collisions, and made a very precise measurement of the Rb-^He enhancement factor. Future use 
of ^^Ne in precision measurement applications would benefit from the use of contact shift polarimetry, but there is 
at present no published measurement of the Rb-^^Ne enhancement factor. In this paper we present the first such 
measurement. We measured the net Larmor frequency shifts, simultaneously impressed by the presence of polarized 
Rb, on ensembles of polarized ^He and ^^Ne precessing in the same [nearly] spherical glass cell. The measurement 
presented in this paper is not the first in which a ratio of contact shift enhancement factors was determined. Baranga 
et. al. measured the shifts of the Larmor precession frequencies of Rb and K vapors residing in the same cell, induced 
simultaneously by the presence of polarized ^He jlH . 

The experiment is described in Section II. Determining the ratio of the Rb-induced shifts permitted us to divide 
out the dependence on Rb magnetization, and express the ■^^Ne enhancement factor K21 in terms of the '^He enhance- 
ment factor K3. The very accurate measurement of K3 by Romalis and Gates determined the absolute contact shift 
enhancement factor K21. Analysis of the data is presented in Section III. Suggestions for improved measurements are 
given in Section IV. 

II. APPARATUS, PROCEDURE, SYSTEM PARAMETERS 

We employed a ^^Ne/^He dual noble gas Zeeman maser to establish simultaneous precession of the ^^Ne and the ■^He 
in the same sealed glass cell. The ^^Ne/'^He maser is a new device that will be described elsewhere; that description 
is summarized here. Figure |l| is a schematic representation of the maser. Gircularly polarized light resonant with the 
Rb Dl transition (795 nm) is generated by a 30 W laser diode array (LDA); the light is directed onto a 1.4 cm dia 
[near-] spherical sealed glass cell containing 228 torr of ^^Ne, 266 torr of ^He, 31 torr of N2, and Rb metal of natural 
isotopic composition |Q. The LDA light propagates along the axis of a uniform magnetic field of strength 3 Gauss. 
The field is generated by a precision solenoid//^- metal-shield system capable of producing a 3 G magnetic field with 
< 35 /iG/cm magnetic field gradients, with no gradient shimming. The main solenoid was driven using a current 
source capable of ~ 3 ppm stability over timescales of order 1000 sec. The glass cell was heated both by the LDA 
light and blown hot air to a temperature of (124.8 ± 1.9)°G (see Sec. III). The temperature of the air surrounding 
the cell was controlled to about 43 mK RMS. The total output power of the LDA was actively controlled to about a 
part in 10^. 

The Rb was polarized by optical pumping; the resulting Rb electron polarization was very close to unity. The output 
power of the LDA was about 30W, so that even though it had a very broad (~ 1.5 nm) linewidth, the on-resonant 
optical power was roughly 0.5 W. Thus, the optical pumping rate was > 10^ photon absorptions per sec per Rb atom. 
The rate of Rb polarization destruction events per Rb atom in the cell was < 10^ sec"^ so that the Rb electron 
polarization was greater than 0.99. Pickup coils wound in series were placed in close proximity to the cell. The coils 
were excited by the magnetic fiux associated with the precession of the polarized noble gases. The coils were part of a 
tuned resonant circuit with resonances at both noble gas precession frequencies. The large current fiow induced in the 
pickup coils by resonant excitation resulted in magnetic fields oscillating at the noble gas nuclear Larmor frequencies 
being impressed back on the precessing atoms. This feedback resulted in steady state maser operation. 

^^Ne has nuclear spin 3/2, and thus a quadrupole moment. The quadrupole moments of the ^^Ne ensemble interact 
coherently with electric field gradients at the glass cell walls, shifting the four Zeeman energy levels so as to split 
the three [otherwise degenerate] dipole transitions |2^. The splitting is proportional to 5ujq oc 3cos^6' — 1, where 
9 is the angle of the cell 's sy mmetry axis with the magnetic field. We oriented the cell as close as possible to the 
"magic angle" cos 9 = •\/l/3 to minimize the quadrupole splitting. The quadrupole wall shifts in the cell were not 
measured; however, computer simulations indicate that without proper cell orientation, realistic-sized splittings could 
disrupt ^^Ne maser oscillation. In any case, ^^Ne maser oscillation with steady-state amplitude was attained. We 
observed a small (~ 10 /iHz) sinusoidal modulation of the ^^Ne maser frequency. The period of the modulation was 
about 16 hr. We speculate that the modulation was a consequence of remnant quadrupole splitting resulting from 
cell orientation error. In any case, the effect was small, with a period much longer than the duration of our contact 
shift measurements, and thus was not a source of significant systematic error. 

The '^He and ^^Ne maser signals were presented to a low-noise preamplier. The preamp output was then analyzed 
using lockin amplifiers, which were referenced to a set of synchronized signal generators tuned near the noble gas 
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maser frequencies and locked to an H-maser-derived 5 MHz signal. The phases of both noble gas masers' precession 
were recorded digitally at a sample rate of 1 Hz. The A/D sampling trigger was also derived from the master 5 MHz 
reference. 

The measurements consisted of observing the maser phases during a series of reversals of the direction of the Rb 
polarization. The Rb polarization was reversed by rotation of the quarter-wave plates (Fig. 1) and was typically 
achieved in < 20 sec. The total duration of a set of reversals was about 1000 sec, limited by magnetic field stability 
as well as the extent to which the noble gas masers' steady states were disrupted by the alternating Rb polarization. 
Two sets of reversal data were acquired. Clear changes in the slopes of the phase curves were precisely correlated with 
reversal of the Rb polarization direction; these frequency changes were the result of the contact shift combined with 
the [classical] shift due to the Rb magnetization distribution. The change of noble gas maser oscillation frequency 
induced by Rb polarization reversal was about 44 mHz for ^^Ne maser and 71 mHz for the ^He maser. 



III. DATA AND ANALYSIS 



Phase data for the second of the two measurement scans are shown in Figs. |^ and||. We model the ^-'^Ne phase (^21 
(see Fig. ||) and the '^He phase (^3 (see Fig. ^ acquired by the lockin amplifiers as evolving according to the following 
equations: 

<f21 = 'P21,o + 721^0^ - 2TTVref,21 + 721 / dt'6B {t') + {k21 - 1 + JdOf) 721 -5" / dt' Mfib {t') 

Jo «J Jo 

^3 ^ y^3,o + IsBot - 2T:Vref.3 + 13 dt' 6 B {t') + {k^ - I + f dDf) l3^ dt'MRb{t') (4) 

Jo «J Jo 

where: 721 (73) is the gyromagnetic ratio of Ne (He); K21 (^3) is the enhancement factor for Ne (He); Jddf is the 
dimensionless factor to account for the classical magnetic field due to the Rb magnetization; ^'re/,2i and Vref,3 are the 
reference frequencies presented to the lockin amplifiers; Bq is the time-average magnetic field strength over the time 
interval; SB is the time-dependent deviation of the magnetic field from its mean value; and Mf/h is the Rb magneti- 
zation. Note that independent maser frequency measurements determined the ratio of the two gyromagnetic ratios, 
73/721, to an accuracy much better than 1 ppm. Also, the effects of magnetic fields due to noble gas magnetization 
fields on the phases have been neglected in comparison to the contact shift effects; this is reasonable for the noble gas 
magnetizations of this experiment (see below). 

Magnetic field drift occurred over the course of the measurement scan. However, the ^^Ne and "^He masers are 
co-located magnetometers, so that we can construct a data set from which the effects of magnetic field drift and Rb 
distant dipole fields have been removed. We define 

Aip = (f3- —(p2i (5) 
721 

It is easily shown that 

/ 73 \ Svr /■* 

Aif = ips^o - V21,o + 27r l're/,21 - '^ref.3 ] t + {K3 - K21) l3^ I dt'Mnb {t') (6) 

V721 / i Jo 

Aside from a trivial phase offset and known linear phase evolution, A(/3 refiects the action of the Rb contact shifts only: 
phase evolution due to magnetic fields is subtracted off. This combined phase profile is plotted in Fig. ^. As seen in 
Fig. ^, the phase evolution due to the contact shift was a piecewise- linear "sawtooth" profile. Since our measurement 
is derived from computing the differences in slopes between adjacent linear regions of the curve, removal of constant 
and linear dependence from the phase data had no effect on the measurement. Thus, our plots of the phase data 
shown in Figs. |^, and ^ are displayed with mean value and mean slope of zero. 

If the Rb magnetization were known, one could immediately deduce K21 using eqn (|^). Even though the Rb electron 
spin polarization was very close to unity, the density was not known a priori. It was therefore necessary to consider 
the ratio of contact shifts impressed on the ip2i and the A(p phase curves, and then determine the temperature 
self-consistently. 

The data were analyzed in two ways which yielded essentially identical results. First, we computed the average 
changes in slope of the phase curves ip2i and Ayi, caused by the reversals of the Rb magnetization. That is, a least- 
squares linear fitting routine was used to extract the slope of each piecewise linear region of the phase curves. The 
differences in slope of adjacent regions was then recorded for each of the phase curves. The mean absolute values 
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of the slope differences were then computed for each data set. We define M21 as the mean of the absolute value of 
the slope differences for ip2i, and AM as that for the Acp profile. The uncertainties were estimated as the standard 
deviation of the mean in the slope differences' absolute values. 

We define the ratio of mean slope changes as TZ = M21/AM. The weighted average value for TZ from the two 
data sets was TZ ~ 0.1262 ± 0.0025 where a one-sigma uncertainty is reported. The uncertainty was almost entirely 
due to that in M2i- The fractional uncertainty in AM for both of the data sets was < 0.1% This indicates that Rb 
polarization reversal was achieved to high precision and that the Rb magnetization was stable to better than 0.1% 
over the course of the measurements. The fractional uncertainty in M21 was much larger, 2%, caused by the effects 
of magnetic field drift. In principle, we could have analyzed the tpa profile rather than (^21, but magnetic field drift 
effects were 73/721 = 9.65 times larger on (^3 than ip2i, making it impractical to extract slope data from ip^ directly 
(see Fig. |). 

In our analysis we did not account for possible variations in frequency due to changes in the noble gas longitudinal 
polarization that occurred in the course of a measurement. The magnetization of the noble gases created magnetic 
fields which varied in time during the course of the each of the two measurement scans, as evidenced by [small] 
maser amplitude changes observed during the scans. A given maser's magnetization field could not affect its own 
frequency, since magnetization fields were parallel to the magnetization in the near-spherical cell. However, each 
maser's frequency was shifted by the other's magnetization, and such shifts did not subtract out of Aip as did 
magnetic field drift. The noble gas polarizations underwent both slow drift as well as variations correlated with 
the Rb polarization reversals. The small uncertainty in the Aip slope differences proves that frequency variations 
due to slow drift in the noble gas magnetizations had negligible effect. We also determined that the variations of 
the noble gas polarization caused by the periodic Rb polarization reversal had negligible effect on the contact shift 
measurement. Periodic reversal of the Rb polarization (i.e., "square wave" modulation of the Rb polarization) induced 
a "sawtooth wave" modulation (i.e., integrated square wave) on the noble gas polarizations of peak-to-peak amplitude 
approximately jset, where ^se is the Rb-noble gas spin exchange rate per noble gas atom, and r = time between 
Rb polarization reversals. Thus, the varying noble gas polarizations induced large changes in the maser frequencies 
over a measurement period: we estimate that a 2.8 mHz shift was induced on the ^He maser by the ^^Ne polarization 
variation and a 0.7 mHz shift was induced on the ^^Ne maser by the '^He polarization variation. However, given a 
constant time interval r between Rb polarization reversals, the average frequency shifts due to the Rb-driven noble 
gas polarization variations did not change from one period to the next, thus they subtracted out in the calculation 
of the period-to-period frequency changes M21 and AM . We conclude that noble gas polarization- induced magnetic 
fields had no systematic effect on the measurement of 7?,. 

It is easy to show from eqns M) and (o) that the ratio TZ can be interpreted as 
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721 K2I + if DDF - 1) 



73 H3 - K21 



(7) 



where it is seen that dependence on the Rb magnetization divided out. The value and estimated uncertainty for TZ 
can then be used to compute the value and uncertainty of the desired contact shift enhancement factor for ^^Ne, 
K21, though the system temperature must be known because ^He contact shift K3 is temperature-dependent. K21 and 
the system temperature must be determined self-consistently, and the effect of classical magnetization fields must be 
accounted for: these questions are addressed below. 

We now consider the second method of analyzing the data to obtain a value for TZ, which yielded essentially the 
same results as the above approach. While this second method is slightly more complex, it is useful because it 
provides a framework for estimating the size of magnetic field drift in the present measurement. This second method 
of extracting a value of TZ from the data considers the phase data sets to be iV-component vectors (denoted by a 
tilde), where N = number of data in the set (~ lO'^ for each of the two scans). We use the Aip profile as a model 
of the shape of the Rb magnetization-induced phase (see eqn ^). We then seek to describe the <^2i profile vector 
as a linear combination of elements of a partial basis of orthonormal vectors. The partial basis included an element 

corresponding to the [normalized] A^p profile (= Vo). Additional basis vectors orthogonal to it, \ Vi> were 

derived from Aip and Legendre polynomials {^'i}i=o at^-i using Gram-Schmidt orthogonalization The number 
of basis vectors Ni, was of order 10^, whereas the dimension of the vector space TV was of order 10'^, so Ni, << N, and 
we find: 

'f2i = ^ aiVi (8) 

i=0 

The coefficient of Vq in this linear combination, oq, was proportional to the value for TZ obtained from that data set. 
The coefficients of the linear combination were determined using linear least squares fitting and by simple projection 



4 



of (pi\ onto the vectors of the partial basis, with identical results. The number of elements in the basis set was varied, 
but the coefficients obtained from the least squares fit remained the same regardless of the size of the basis set. 

Using this second data analysis method, we thus determined two values of 7?., one from each data set. It was 
not possible to estimate uncertainties for the individual values of from the vector analysis method. Uncertainty 
estimates from the least squares analyses were unrealistically small because the vector analysis could not account 
for the effect of a non-zero projection of the magnetic field drifts along the Rb magnetization-induced phase profile 
vectors. We take the value of to be the mean of the results for the two data sets, and we coarsely estimate the 
uncertainty of 7?, as the deviation of the two values, to obtain 7?, = 0.1287 ± 0.0036. This result agrees well with that 
of the first method of data analysis. For computing K21 we will use the result from the first method of analysis since 
the error estimate is probably more reliable. 

We can compute the magnetic field drift profile by subtracting from (p2\ its projection on the Rb magnetization 
basis profile Vq: 

= :^^ (^21 - (<^2i -14)14) (9) 

Note that the profile 8B is the magnetic field drift orthogonal to Vq. The drift profile from the second data set is 
plotted in Fig. ||; the time derivative of the [discretely sampled] profile was estimated by finite difference. Note that 
the field was stable to about ±3 ppm RMS. 

Having obtained a value for TZ, we can then solve eqn (|^) for K21 in terms of measured parameters, since K3 is known 
However, K21 is expected to be temperature-dependent [|l^,0 and so the temperature at which the experiment 
was carried out must be reported along with the value of K2i. The Rb and noble gas temperature indicated by the 
hot air temperature control system was in error, because: the absolute calibration was poorly known; and the cell 
was heated by the LDA to a temperature higher than that indicated by the control system. Thus, it was necessary 
to deduce the noble gas/Rb vapor temperature from the contact shift data, and the known "^He enhancement factor 
(with its known temperature dependence), in a self-consistent way. This was achieved by first noting that the mean 
slope change AM is related to the enhancement factors and the Rb magnetization via 

AM = 273y |(k3 - K21) Mr, (T)| (10) 

where {Mm, (T)\ (see eqn (||) above) is a known function of temperature. The Rb number density [Rb] is related to 
temperature by |26[ 

-|^q9.318-4M0 

[^^]= (1.38.10-^V).^ c..-3 iT^n^K) (11) 

The measured temperature dependence of the -^He enhancement factor is 

K3 = 4.52 + 0.00934 [r(°C)] (12) 

Substituting ( |lT|) and (|^) into (10) yields an equation relating K21 and temperature (and known/measured parameters). 
Substituting ^^) into (0) yields another equation relating K21 and temperature. This system of two [nonlinear] 
equations and two unknowns can be solved to yield the temperature and K21. The temperature determination 
depends strongly on the enhancement factor, while the enhancement factor determination depends only weakly on the 
temperature (via the known temperature dependence of K3). Thus, an iterative method to compute the temperature 
rapidly converges: an initial guess for the temperature enables calculation of K21, which in turn permits calculation of 
the Rb magnetization and thus the system temperature. This corrected temperature is used to re-compute a better 
estimate for K21, which is used to re-compute a better estimate of the Rb magnetization, etc. 

Using the result TZ — 0.1262 ±0.0025 from the slope difference analysis, we determine the contact shift enhancement 
factor for Rb and ^^Ne to be K21 — 31.8 ± 2.8 at a temperature of (124.8 ± 1.9)°C. Note that the noble gas/Rb vapor 
temperature is much higher than the value of 112°C indicated by the hot air temperature control system, presumably 
the result of heating by the laser diode array. We are not aware of a previous measurement of the Rb-^^Ne contact shift. 
Walker estimated the Rb-^^Ne enhancement factor to be 38, and the ratio of the Rb-^^Ne and Rb-^He enhancement 
factors to be 4.3 ]Q. We report this ratio to be 5.6 ± 0.5. 

The 3% uncertainty in TZ dominates other sources of error, inducing a 9% uncertainty in K21. The contribution to 
the uncertainty in K21 resulting from the 1.5% uncertainty in the He enhancement factor K3 is only 1.5%, which is 
negligible. The uncertainty in temperature induced a 0.3% uncertainty in K21. To the extent that the cell could be 
described as an ellipsoid of revolution, its deviation from sphericality induced an uncertainty of 0.15% in K21 with 
zero systematic correction. 
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The claim of zero systematic correction of the classical magnetization field shift for a ellipsoidal cell is justified 
by noting that: i) the classical magnetization field inside an ellipsoid of uniform magnetization is uniform, with 
its magnitude depending on the ellipsoid's shape and the relative orientation of the ellipsoid major axis and the 
magnetization vector; and ii) when the ellipsoid major axis is oriented at the magic angle 9m — cos" -1 (1/V3) with 
respect to the magnetization direction, the component of the classical magnetization field parallel to the magnetization 
is exactly that of a spherical distribution with the same magnetization, regardless of the ellipsoid's shape. We oriented 
the cell at 9m in order to suppress quadrupole wall shifts, with an uncertainty of about 10°. The 10° uncertainty in 
orientation leads to the estimated uncertainty of 0.15% in K21 due to measured asphericality of the cell. 



IV. SUGGESTIONS FOR IMPROVED MEASUREMENTS; CONCLUSIONS 

There are straightforward improvements that could be implemented to reduce the uncertainty in the present mea- 
surement. The experiment was conducted in the course of developing a ^^Ne/^He dual Zeeman maser, so that 
opportunity for taking contact shift data was very limited. A factor of three improvement in the result could easily 
be achieved just by taking ten times more data. Active field stabilization using an independent magnetometer p^ ] 
would essentially eliminate field drift, reducing statistical uncertainty in the measurement to that due to phase noise, 
a factor of > 0.03/0.001 = 30 reduction. Active field stabilization would also permit a reduction in the phase noise 
since the signal acquisition bandwidth could be reduced by at least a factor of 3, so that the statistical uncertainty 
could be reduced by about a factor of 300 compared to the present experiment. This uncertainty would already be 
less than the uncertainty contributed by cell asphericality (0.15%). However, the absolute uncertainty in the Rb-^-'^Ne 
contact shift enhancement factor K21 would then be limited by the uncertainty in that for He, K3, 1.5%, whereas the 
ratio would be known to an uncertainty of only a few tenths of one percent. 

It should be possible to measure ratios of alkali-noble gas enhancement factors to a level well below a part in 
10'^ if one could completely eliminate effects of the classical magnetization field shift. We have already shown that 
one can use differential magnetometry to obtain data from which classical magnetization field shifts are eliminated 
to high precision. It can be shown that using a three-species noble gas maser, one could measure the contact shift 
enhancement factor of a third species if the other species' enhancement factors are known, independent of the classical 
magnetization field shift. Even though it has not yet been done, we believe it is a straightforward matter to implement 
a three-species maser in a single-bulb glass cell. 

A comprehensive set of contact shift data would serve as a probe of noble gas-alkali interaction potentials [ pO[ . 
A program using multi-species noble gas masers to measure contact shift ratios could be implemented using ^^®Xe, 
*^Kr, ^^Ne, and '^He along with K, Cs, and Rb in turn. Na would be a poor candidate for use in a noble gas maser 
because of its relatively low vapor pressure. The remaining stable noble gas isotope, ^^^Xe, cannot be used with the 
other noble gases in a multi-species maser because its nuclear dipole moment has a sign opposite to that of the others 
(making it impossible to achieve a population inversion simultaneously with the other species). However, it may be 
possible to measure the ^'^^Xe precession frequency in free induction decay in the presence of two other noble gas 
masers. There are no other stable noble gas isotopes having non-zero nuclear spin. 

Contact shift ratios would be converted to absolute values for enhancement factors by using the measured values 
for ^He/Rb and '^He/K [|l^. The enhancement factor for ^He/Cs has not yet been measured. 

Measuring contact shift ratios to a few parts in 10^ could proceed as follows: one would first measure the ratio 
of the enhancement factors of *'^Kr and ^^^Xe to each other. The contact shift enhancement factors for *'^Kr and 
^^^Xe are estimated by Walker to be about 230 and 730, respectively [Q. Using an independent magnetometer with 
the colocated Xe-Kr masers, the ratio of the [large] enhancement factors of Xe and Kr could be measured to a few 
parts in 10"* even if the classical magnetization field were accounted for only to the 10% level (which can be done 
trivially- a ~ 1% accounting is claimed in this work). Triple masers using Kr and Xe, along with ^^Ne and ^He in 
turn would yield high-precision measurements of the ratios of all the enhancement factors; the known enhancement 
factor for He would then determine the absolute values of the other enhancement factors. Temperature dependences 
could of course also be measured. In particular, enhancement factor ratios as a function of density (normalized to 
some experimental reference density) could be determined to very high precision. Improved cell thermometry would 
eliminate the need to depend on a temperature-density calibration p6| to determine the cell temperature. Also, for 
very high precision measurements, control of the noble gas magnetization fields would have to be maintained more 
carefully than in the present work; in particular, the time between polarization reversals could be decreased in order 
to increase the relative size of the contact shifts in comparison to noble gas polarization induced frequency shifts. It 
might also be important to account carefully for the contribution of van der Waals molecule formation to the contact 
frequency shifts. 

Finally, we remark that a two species maser using ^He and ^^Ne has been proposed for a high precision test of 
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local Lorentz invariance (LLI) in which Rb-noble gas contact shifts would be an important systematic effect. 
Adding a co-located ^^^Xe maser would provide a high precision in situ measurement of the Rb magnetization. With 
part-in-10^ knowledge of the enhancement factors, the systematic effects of contact shifts on a search for LLI violation 
could be rendered negligible with the use of the -'^^^Xe maser " co-polarimeter" . 

We thank T. Chupp for the use of the ^^Ne/'^He cell, and gratefully acknowledge D. Phillips and D. Bear for useful 
conversations. 
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FIG. 1. The ^^Ne/^He dual Zeeman maser. The device operates similarly to the ^^^Xe/'^He maser except that the 
^^Ne/^He maser employs a single- rather than a double-bulb cell. Co-located ensembles of noble gas atoms are polarized by 
spin exchange collisions with optically pumped Rb atoms. Precession of polarized atoms excites current in nearby inductive 
pickup coils, which are part of a circuit tuned with resonances near the Larmor precession frequencies of the two noble gas 
species. The current induces magnetic fields which act back on the precessing atoms, providing feedback for maser oscillation. 
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FIG. 2. Example ^^Ne maser phase profile from one (of two) measurement scans. Phase dependence linear in time has been 
removed. 
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FIG. 3. Example ^He maser phase profile measured simultaneously with the data of Fig. 2. Phase dependences linear in 
time have been removed. Note the considerable distortion in this phase curve due to magnetic field drift. 
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FIG. 4. The combined phase of the two masers as per eqn (^). This profile is proportional to the phase variation induced 
by the Rb-noble gas contact shift alone. The sawtooth shape arises from periodic reversal of the Rb polarization, which reverses 
the Rb-noble gas contact frequency shift. The highly regular shape of the difference curve reflects the fact that magnetic field 
drift effects are subtracted out. 
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FIG. 5. Magnetic field drift during the measurements of Figs, g and q, computed from eqn The field strength was 3 G, 
so that the fractional RMS stability of the magnetic field was about 3 ppm. Field drift was the dominant source of experimental 
error. 
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